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ABSTRACT 



Set of neutron star observational results is used to test some selected equations of state of dense 
nuclear matter. The first observational result comes from the mass-baryon number relation for pulsar 
B of the double pulsar system J 0737-3039. The second one is based on the mass-radius relation 
coming from observation of the thermal radiation of the neutron star RX J 1856.35-3754. The third 
rG • one follows the population analysis of isolated neutron star thermal radiation sources. The last one 

Mh' is the test of maximum mass. The equation of state of asymmetric nuclear matter is given by the 

^ , parameterized form of the relativistic Brueckner-Hartree-Fock mean field, and we test selected pa- 

$_( ■ rameterizations that represent fits of full relativistic mean field calculation. We show that only one 

^ ' of them is capable to pass the observational tests. This equation of state represents the first equation 

C^ , of state that is able to explain all the mentioned observational tests, especially the very accurate test 

given by the double pulsar even if no mass loss is assumed. 
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Tj- . 1. Introduction 



Neutron stars are compact objects that play important role in different areas 

of modem physics. Here we concentrate our attention on the possibility that phe- 

^ I nomena related to neutron stars can be used as tests of equation of state (EoS) of 

asymmetric nuclear matter. The tests used in this paper represent a subset of tests 

used previously by Klahn et al. (2006). We focus our attention on tests that come 

k><( ; from astronomical observations. However we have not applied the very promising 

^ ■ test coming from the observations of quasiperiodic oscillations (QPOs), since the 



theory and data interpretation is still in progress (see e.g. Torok et al. (2008a,b,c, 
2010) van der Klis (2004)). The QPO test apphed on the 4U 1636-536 object in 
Klahn et al. (2006) represents the maximum mass test in the present paper, and 
another test following the QPO phenomena observed in the 4U 0614+09 object do 
not provide a strong test and all the EoS tested in this paper pass it. 

A wide spectrum of different equations of state of nuclear matter and their 
applications to astrophysical problems has been reported in literature (see, e.g., 
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Haensel, Zdunik and Douchin 2002, Rikovska Stone et al. 2003, Weber, Negreiros 
and Rosenfeld 2007, Lattimer and Prakash 2007, Burgio 2008). Some of the EoS 
collections (even though not all of them are up-to-date already) give an amazingly 
rich general overview of the state-of-the-art, whereas the others emphasize some 
specific aims. All these EoS yield (nearly) the same properties close to the standard 
nuclear density (Pa^ « 0.16 nucleon/fm^ ^ 2.7 x 10^"* g/cm-'), but when one is 
far off this value, s/he has to rely more on underlying principles than on possible 
experimental verification of predicted physical observables. 

Here we concentrate our attention on relativistic asymmetric nuclear matter 
where the EoS stem from an assumed form of the interaction Lagrangian. The 
calculations use the relativistic mean-field theory with allowance for an isospin de- 
gree of freedom (Kubis and Kutchera 1997, MUler and Serot 1996). We employed 
the Dirac-Brueckner-Hartree-Fock mean-field approach in its parameterized form 
suggested in Gmuca (1991) which reproduces the nuclear matter results of Huber, 
Weber and Weigel (1995). That has been used to calculate high-density behavior 
of asymmetric nuclear matter with varying neutron-to-proton ratio (Gmuca 1992). 
The proton fraction has been determined from the condition of P -equilibrium and 
charge neutrality, and it is density-dependent. We have extended our calculations 
for densities up to 4 x p^r and if there was an astrophysical motivation even higher. 

The EoS is used to model the static, spherically symmetric neutron star in the 
framework of general relativity. The equation of hydrostatic equilibrium is solved 
for different central parameters (pressure, energy density, baryon number density). 
The radius of the neutron star model is then given by the condition of vanishing 
pressure. The resulting properties of the neutron star model are then compared with 
observational data. From the test ensemble presented by Klahn et al. (2006) we 
choose four astrophysical observations to test our selected parameterizations that 
have been found to be a good description of nuclear matter at subnuclear densities 
for pure neutron matter and up to 2 x p/v for symmetric nuclear matter (Kotulic 
Bunta and Gmuca 2003). 

The maximum mass test is the standard way to test the EoS of asymmet- 
ric nuclear matter (see e.g. Haensel, Potekhin and Yakovlev 2007, Lattimer and 
Prakash 2007, Klahn et al. 2006). The usual value to constrain the maximal mass 
of neutron star comes from observations of double pulsar PSR 0751-1-1807 giving 
M = (2.1 ±0.2) Mq, with Mq being the solar mass. This value was, however, 
lowered to M = (1.26±0.14) Mq (Nice, Stairs and Kasian 2008) and could not be 
used as maximum mass test anymore. Another value that could serve as maximum 
mass test comes from the observation of QPOs. The mass is constrained on the 
basis that the observed frequency corresponds to the frequency at innermost stable 
circular orbit (Barret, Olive and Miller 2005, Belloni, Mendez and Homman 2007, 
van der Klis 2004). 

Popov et al. (2006) used the population synthesis of the isolated neutron star 
sources of thermal radiation and concluded that the neutron stars with mass M < 
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1 .5 Mq could not cool via direct URCA reactions. This conclusion follow from the 
fact that all observed sources of thermal radiation have masses bellow the quoted 
value. This could be explained by the fact that more massive objects cool via 
the direct URCA reactions which represents the fast cooling scenario and thus the 
thermal radiation could not be detected. These arguments were used by Klahn et 
al. (2006) to build a strong and a weak test on EoS. 

A very accurate test of EoS was developed by Podsiadlowski et al. (2005). They 
based it on the model of double pulsar system J 0737-3039 formation. The model 
predicts the pulsar B of this system to be born via the electron capture supernova 
what suggests extremely low mass loss and thus the number of particles conserved 
during the progenitor collapse to neutron star. This put limits on the mass-baryon 
number relation. Instead of the baryon number that represents the total number of 
baryons contained in the neutron star, the baryon mass could be used equally. 

The thermal radiation coming from the neutron star source RX J 1856.35-3754 
could be used to put limits on the mass-radius relation of the neutron star model. 
TriJmper et al. (2004) used two different models to explain the spectral feature for 
this specific source and found its apparent radius that represents the radius of the 
neutron star as seen by a distant observer. The analyses of data to obtain the isolated 
neutron star radius strongly depend on the radiation spectrum emitted by the object 
and the estimated radius is proportional to the distance from Earth to the source. 
The distances obtained for RX J1856.5-3754 range from D = 6llg pc (Walter and 
Matthews 1997) to D = 161+}4 pc (van Kerkwijk and Kaplan 2007). The derived 
apparent radius /?oo is given by the model of the atmosphere. The original model 
by Pons et al. (2002) resulted in R^/D = 0.13 km.pc^^ . Trumper et al. (2004) 
presented new models of atmosphere leading to the estimates of /?co = 16.5 km 
for the two component model of spectra and R^ = 16.8 km assuming continuous 
temperature distribution model. If the distance derived by van Kerkwijk and Kaplan 
(2007) and the original model of Pons et al. (2002) are used together, they lead to 
unexpectedly high estimate /?oo = 20.9 km. Recently Steiner, Lattimer and Brown 
(2010) presented results based on new analysis of data giving the distance 119 ± 
5 pc and the original model for atmosphere (Pons et al. 2002) then implies /?„ = 
15.47 km. We decided to use the three values Roo = 15.5, 16.8, 20.9 km to put 
limits on neutron star equation of state. 

Another promising way to constrain the equation of state are the moment of 
inertia measurements (see e.g. Lattimer and Prakash 2007 and references therein). 
Two ways have been proposed quite recently. One for the Crab pulsar (Bejger 
and Haensel 2002,2003) following observations of the pulsar-nebula system, and 
the other for the pulsar A of the double pulsar system J0737-3039 (Bejger, Bulik 
and Haensel 2005) based on the measurements of the second order post Newtonian 
parameters of the binary system. Even thought both ways could provide strong 
limits on the equations of state in principle, they need more accurate observational 
inputs. We need better estimates of the mass of Crab nebula in the first case and 
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very accurate measurements of orbital parameters are necessary to calculate the 
moment of inertia in the second case. For these reasons we do not include these 
tests to our calculations. The measurements of moment of inertia of the neutron 
star together with its mass put limits on the radius of the neutron star that is crucial 
for the cooling scenarios (see e.g. Lattimer and Prakash 2007, Stuchlik et al. 2009). 
The paper is organized as follows. In section 2 we present our EoS and details 
of the neutron star matter description. Section 3 briefly summarizes the model of 
static spherically symmetric neutron star. We present our results and compare them 
to observations in section 4. The paper is closed by conclusions in sections. 

2. Equation of state of neutron star matter 

2. 1. Asymmetric nuclear matter in relativistic mean-field approach 

We follow the Dirac-Brueckner-Hartree-Fock (DBHF) mean field (see Weber 
1999, Walecka 2004, de Jong and Lenske 1998, Krastev and Sammarruca 2006 
for underlying theories), which easily allows to consider different neutron-proton 
composition of the neutron star matter, and also the inclusion of non-nucleonic 
degrees of freedom. 

The full mean-field DBHF calculations of nuclear matter (Huber, Weber and 
Weigel 1995, Lee et al. 1998, Li, Machleidt, and Brockmann 1992) have been pa- 
rameterized by Kotulic Bunta and Gmuce (2003), and we employ their parameter- 
ization with one-boson-exchange (OBE) potential A of Brockmann and Machleidt 
Li, Machleidt, and Brockmann (1992). We refer to the paper of Kotulic Bunta 
and Gmuca (2003) for the explicite set of values of the corresponding parameters. 
The model Lagrangian density includes the nucleon field \|/, isoscalar scalar me- 
son field a, isoscalar vector meson field (O, isovector vector meson field p ,and 
isovector scalar meson field 5, including also the vector cross-interaction. The 
Lagrangian density in the form used by Kotulic Bunta and Gmuca (2003) reads 

x (\|/, a, CO, p , 5) = \i/[Y^(/a'' - gcoo/') - (mN - gca)]\|/ 
+-(a^aa^a - m^^o^) - -co^vo/"' + -m^^fn^af 



-T^aWTNfeaCT) - -Ca{ga<^) + ^^wfem CO^CO") 
+-(3^53^5 -m8252) + -mp2p^p^--p^^p 



.A-V 



+:^Av(gp2p^p'')(gffl^(0^a/') -,?pp^?'/'x\^ + 585\^x\^, (1) 



2 
where the antisymmetric tensors are 



CO^v = 3vC0^ — 3^C0v, 

P^v = 3vP^-3a,Pv; (2) 
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the strength of the interactions of isoscalar and isovector mesons with nucleons is 
given by (dimensionless) coupling constants g's and the self-coupling constants 
(also dimensionless) are b^ (cubic), c^ (quartic scalar) and Cq, (quartic vector). 
The second and the fourth lines represent non-interacting Hamiltonian for all me- 
sons, Ay is the cross-coupling constant of the interaction between co and p mesons. 
Furthermore, wtn is the nucleon mass, 9^ = g^ and y's are the Dirac matrices 
(Kotulic Bunta and Gmuca 2003, Serot and Walecka 1986, Weber 1999). 

We choose here three following parameterizations, which were shown to yield 
the best fits to the well-known properties of nuclear matter 

H HA in Kotulic Bunta and Gmuca (2003) represents the best RMF fit to results 
obtained by Huber, Weber and Weigel 1995. 

L LA in Kotulic Bunta and Gmuca (2003) represents the best RMF fit to results 
obtained by Lee et al. 1998, but does not include the 8 mesons to nucleons 
coupling. 

M MA in Kotulic Bunta and Gmuca (2003) represents the best RMF fit to results 
obtained by Li, Machleidt, and Brockmann 1992, but does not include the 5 
mesons to nucleons coupling. 

The EoS of Kotulic Bunta and Gmuca which have been found to be a good 
description of asymmetric nuclear matter, are easily expressed up to about 4 x p/v 
(parameterization //) or even higher (parameterizations L and M). 

2.2. ^-equilibrium 

The total energy density of n-p-e-ju matter is given by 

'E='EB{nB,Xp) + 'Ee{ne) + 'E^{n^), (3) 

where Eb('^B5-^p) is the binding energy density of asymmetric nuclear matter, «, is 
the number density of different particles (/ = n, p, e, /n), «b = «p + «n is the baryon 
number density and Xp = «p/«B is the proton fraction. The leptonic contributions 
"Eiini) {I = e, /^) to the total energy density are given by 

Pf{i) 
■Li{ni) = -^ j {my + p^cY^4np^dp, (4) 



where pf^^ is the Fermi momentum of Z-th kind of particle. 

The matter in neutron stars is in P -equilibrium, i.e. in equilibrium with respect 
to nop + e^ op + /i. The (anti)neutrinos contribution could be neglected, be- 
cause the matter is assumed to be cold enough that they can freely escape. The 
equilibrium is given by equality of chemical potentials /Ha = /^p + /7e = /7p + /^^ , 
where the chemical potential of each kind of particle is given by /i, = d'E /drii . 
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The chemical potentials of electrons and muons are simply fA\ = ^/rn^c'^ + pj^, 
while the chemical potentials of nucleons are 



^r 



""(p. n) 

The binding energy density of asymmetric nuclear matter could be expressed 
in terms of proton fraction Xp (Danielewicz and Lee 2009) 

£B(nB,-^p) = £snm(«b) + (1-2xp)^5'(?ib), (6) 

where Esnm is the energy density of symmetric nuclear matter (xp = 0.5) and 
5'(«b) is the symmetry energy density, that corresponds to the difference of binding 
energy density between pure nuclear matter and symmetric nuclear matter 

The symmetry energy 5'(«b) is the factor corresponding to the second order 
term in expansion of binding energy density in terms of asymmetry parameter 8 = 
(?in — «p)/(«n + «p) = 1 — 2xp and reads 



^^^^^_ia^£B(«B,5) 



2 352 



(7) 

S=() 



From equation (6) one can see that symmetry energy is the difference of binding 
energy per particle between pure nuclear matter and symmetric nuclear matter. 

5(?ib) = Eb(«b,-^p = 0) - •Eb(«b,-^p = 0.5). (8) 

The condition of (3 -equilibrium then reads 

i^e = A/^ = iUn - A/p = 4 ^(l-2Xp). (9) 

'^ «B 

and it is solved together with condition of charge neutrality («p = «e + «;u) to obtain 
the proton fraction of neutron star matter. The binding energy per baryon in depen- 
dence on the baryon number density is illustrated in Figure 1 . The proton fraction 
of matter at the beta-equilibrium is given, for the chosen three EOS parameteriza- 
tions, as a function of the baryon number density depicted in Figure 2. 

2.3. EoSfor low densities 

The nuclear EoS have been the dominant input for the calculations in the high- 
density region, namely p > lO''* g/cm-' . For lower densities, the EoS used are the 
following: 



Feynman-Metropolis-Teller EoS for 7.9 g/cm-' < p < lO'* g/cm^ where mat- 
ter consists of e^ and 26^^' Feynman, Metropolis and Teller (1949); 

Baym-Pethick-Sutherland EoS for lO'* g/cm^ < p < 4.3 x 10" g/cm^ with 
Coulomb lattice energy corrections Baym, Pethick, and Sutherland (1971); 
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Figure 1: Binding energy per particle of different types of nuclear matter for used 
parameterizations. Left Matter at P -equilibrium, Middle symmetric nuclear matter, 
and Right pure neutron matter 

• Baym-Bethe-Pethick EoS for 4.3 g/cm^ x 10^^ < p < 10^"* g/cm^: here, 
e ^ , neutrons and equilibrated nuclei calculated using the compressible liquid 
drop model Baym, Bethe, and Pethick (1971). 



3. Neutron star models 



We consider static spherically symmetric models of neutron stars. The inte- 
rior spacetime is described by the internal Schwarzschild metric (see, e.g., Misner, 
Thorne and Wheeler 1973, Haensel, Potekhin and Yakovlev 2007) that can be writ- 
ten in geometrical units (c = G=l)as 



ds^ = -e^^dt^ + e^^dr^ + r^ (de^ + sin^ ed(|)2) , 



(10) 



where the radial component of metric can be expressed as a function of energy 
density p 

„2X_ r 



-, m{r) =471 / prjdri 
Jo 



(11) 



r — 2m(r) 

The matter is assumed to be perfect fluid described by the energy momentum 
tensor 

T^"' = {P + pXu'' + Pg^\ (12) 
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Figure 2: Proton fraction of matter being at (3 - equilibrium for used parameteriza- 
tions. Also lines of direct URCA threshold (marked with xdu) for all parameteri- 
zations are depicted. 

(Misner, Thorne, and Wheeler 1973). where P is the pressure, u^ is the 4 -velocity 
of matter and g^^ is the metric tensor. The energy momentum tensor satisfies the 
conservation law T^^-y = 0. 

The hydrostatic equilibrium is in general relativity given by the Tolman-Oppen- 
heimer-Volkoff equation (TOV) (Oppenheimer and Volkoff 1939, Tolman 1939), 
which reads 

— = -p + P ) ' (13) 

dr r[r — 2m[r)) 

Integration of TOV starting from given central energy density Pc uses the EoS and 
finally yields the radius R, given by the boundary condition P{R) = 0, and the 
gravitational mass M = m{R) of the neutron star. 

Another useful quantity to calculate is the so-called baryonic mass Mb that 
represents the total number of baryons contained in the neutron star multiplied by 
the atomic mass unit u. The baryonic mass is then expressed as 

R 

2m (r) 



Mp 



: 471U / «B [r) 




1 



-1/2 



r^dr, 



(14) 



where nB{r) is the baryon number density at the radius r. 
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Figure 3: Mass given as a function of central baryon number density for different 
parameterizations. The stars correspond to the minimum mass of a neutron star that 
could cool via direct URCA reactions. 

4. Results versus observations 



Several dozens of neutron stars and/or similar objects have their masses re- 
ported; a great majority of them is in very close vicinity of 1.4 Mq, and only 
very few are significantly above (see, e.g., the compilations in Bethe, Brown and 
Lee (2007), Lattimer and Prakash 2007) and observations and analyses (see, e.g, 
Rikovska Stone et al. 2003, Weber, Negreiros and Rosenfeld 2007, Podsiadlowski 
et al. 2005, Triimper et al. 2004, Pons et al. 2002, Kramer and Wex 2009, Krastev 
and Sammarruca 2006, Lattimer and Prakash 2007, Blaschke, Klahn and Sandin 
2008, Dexheimer, Vasconcellos and Bodmann 2008, Klahn et al. 2006, Nice, Stairs 
and Kasian 2008, Rikovska Stone et al. 2007). However, recent results of the data 
fitting of kHz quasiperiodic oscillations observed in the low-mass X-ray systems 
containing neutron stars indicate relatively high masses of M > 2 Mq (Belloni, 
Mendez and Homan 2007, Torok et al. 2008a,b,c, Barret, Olive and Miller 2005, 
Boutelier et al. 2010, Boutloukos et al. 2006) which could provide very strong 
constraint on the EoS. On the other hand, modification of the characteristic orbital 
frequencies by a magnetic repulsion caused by the interaction of slightly charged 
matter in accretion disc in vicinity of a neutron star with dipole magnetic field could 
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shift the mass estimates to lower values close to canonical 1.4 Mq (Bakala et al. 
2008). Our calculations with parameterization H allow for the existence of neutron 
stars even for so heavy masses. 

4.1. Direct URCA constraints 

The proton fraction Xp of matter in (3 equilibrium is presented in Figure 2 to- 
gether with the direct URCA threshold. The direct URCA reactions n ^ p + e^ + 
Vg could operate only if the proton fraction exceeds the threshold given by the 
condition 

■^DU = ; —, (15) 

1+ l+x'J' 



where Xg = ne/{ns + nfj) . One can see that only parameterizations L and M enable 
rapid cooling. The threshold densities are «du = 0.457 fm^'' in the case of pa- 
rameterization L and wdu = 0.571 fm^^ in the case of parameterization M . These 
values correspond (see Figure 3) to neutron star masses M = 1 .47 Mq (parameter- 
ization L) and M = 1.39 Mq (parameterization M). Parameterizations L and M 
thus do not fulfill the direct URCA constraints, however Klahn et al. 2006 used also 
the value 1.35 Mq as a weaker test that is passed also by parameterizations L and 
M. 

4.2. Maximum mass 

The maximum mass limit is probably the most often used test of the equation of 
state. The maximum masses given by EoS used in this paper are M^^^ = 2.18 Mq , 
^max — 1-92 Mq and M^^^ = 1.62 Mq . The maximum mass obtained for objects 
containing matter described by parameterizations L and M follows the require- 
ments of stability with respect to radial oscillations (dM/dnc > 0). In the case H 
we used the values corresponding to the central density «b (r = 0) = 0.66 fm^^ , 
because for the densities above this value the model used for the EoS is not with- 
out questions and also because only with central densities up to about 4 x normal 
nuclear density we were able to explain masses of neutron stars that meet the obser- 
vational requirements. With some extrapolations, higher masses could be in prin- 
ciple modelled, but we decided to use parameterization H up to the quoted density 
only since there is no current astrophysical observation of such a high mass. The 
observation of high mass is however crucial and very promising issue of astrophys- 
ical observations. It should be noted that Klahn et al. 2006 used the value that 
could not be used anymore. Also the result for the source 4U 1636-536 that gives 
M = (1.9-2.1) Mq as proposed by Barret, Olive and Miller (2005) should be used 
rather as an upper limit of the neutron star mass than as its estimate, see, e.g. Miller, 
Lamb and Psaltis (1998) for underlying theories. The neutron star mass is inferred 
due to the highest observed frequency of QPOs observed in the system, under the 
assumption of identifying the highest frequency with the Keplerian frequency of 
the innermost stable circular orbit (ISCO). Clearly this gives an upper limit on the 
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Figure 4: Relation of calculated gravitational mass M and the baryonic one Mb 
for different parameterizations. The limitations imposed by the analysis of the 
J0737-3039 double pulsar are drawn as a small rectangle. The box is also extended 
to the left by O.OOSM© indicating the possible mass loss. 

mass, and the real neutron star mass has to be expected smaller because the QPOs 
have to be excited above ISCO. Up to date, one of the two pulsars Ter 5 I and J has 
a reported mass larger than 1.68 M© to 95% confidence level (see, e.g. Lattimer 
and Prakash (2007) and references therein^. Champion et al. (2008) predicted 
mass of PSR J1903+0327 to be M = 1.74 ±O.O4M0. Freire (2009) estimated the 
mass for the same source to be M = 1.67 ±0.01 Mq. These values, even if they 
are different, give approximately the same limit on mass when they are combined 
together, namely > 1.66 Mq at 2 a level. These predictions are not in favour the 
parameterization M with M^^^ = 1.62 Mq . 



4.3. Double pulsar J07 37-3039 

Podsiadlowski et al. (2005) investigated possible formation scenarios of double 
pulsar J0737-3039. They have shown that one can test EoS assuming the pulsar B 
is formed by an electron-capture supernova. Such scenario enables formation of the 



'The individual pulsar masses unfortunately are not assumption-independent. In our discussion, 
we adhere to the value 1.68 Mq reported by Lattimer and Prakash, but bearing in mind the possible 
uncertainty in its derivation. 
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Figure 5: Mass-radius relation for different parameterizations. The lines corre- 
sponding to RX J 1856.5-3754 gives a lower mass limit, that should the given EoS 
get over. 



pulsar B that has low but very accurately measured mass M = 1 .2489 ± 0.0007 Mq 
(Kramer and Wex 2009). If this pulsar is born under the presented scenario, its 
baryonic mass Mb should be in the range 1.366 to 1.375 Mq. The authors also 
argue the matter loss being low (the matter loss they give is few times 10^^ Mq). 
The relation between the gravitational and the baryonic masses together with the 
limitations derived from the double pulsar observations are presented in Figure 4. 
One can see that the only parameterization that meets requirements assuming no 
mass loss is the parameterization H . The parameterization M is able to explain the 
results if one includes mass loss predicted by Podsiadlowski et al. (2005). Unfor- 
tunately this parameterization was ruled out by the maximum mass test. 

4.4. Isolated neutron star RX J 1856.5-3754 

Several authors (see, e.g., Triimper et al. 2004, Pons et al. 2002, van Kerkwijk 
and Kaplan 2007, Steiner et al. 2010) discussed observations of the isolated neutron 
star RX J 1856.5-3754 and they found constraints on the mass-radius relation of this 
particular neutron star. They found the limits of the apparent radius being given byt 
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the mass-radius relation 

M _ R / R^\ 
M^ ~ 2.95 km V " ^y ' 

that could serve as a test of equation of state. We have used three different val- 
ues for Rco namely Rco = 15.5, 16.8, 20.9 km. None of tested parameterizations 
is able to explain the apparent radius Rc^ = 20.9 km. Parameterization H is the 
only one capable of explaining the apparent radius /?„ = 16.8 km estimated by 
TriJmper et al. (2004). The lowest predicted apparent radius could be modeled by 
all parameterizations considered in this paper. The mass-radius relations for all 
parameterizations together with observational limits are illustrated in Figure 5. 

5. Conclusions 

We have employed the parameterized form of the relativistic mean-field EoS 
for asymmetric nuclear matter with vector cross interaction. The proton fraction 
was varied in accord with the need of the P -equilibrium and charge neutrality. As- 
suming spherically symmetric geometry and using TOV equation, we constructed 
models of neutron stars for different central parameters. We have used set of ob- 
servational data to test EoS of nuclear matter represented by three different param- 
eterizations of relativistic Brueckner-Hartree-Fock equation. 

We have shown that only the parameterization H is able to pass almost all the 
tests considered in this paper. The only exception is the apparent radius /?oo = 
20.9 km estimation for the isolated neutron star RX Jl 856.5-3754; however this 
estimate is based on distance measurements being still widely discussed. This pa- 
rameterization also represents the only EoS based on the relativistic Brueckner- 
Hartree-Fock theory that could explain the formation of pulsar B in the double 
pulsar system J 0737-3039 without mass loss. 

Our present calculations have been done considering only neutrons and protons 
in P -equilibrium with electrons and muons. We aim to continue in tests of given 
EoS in future. One of our plans is to include hyperons. Another is to perform 
more detailed tests based on the promising fitting of observational data of quasi- 
periodic oscillations in low-mass X-ray systems measurements. This necessitates 
to investigate the rotational effects on neutron star models based on the Hartle- 
Thorne metric reflecting mass, spin and the quadrupole moment of the neutron star 
(Hartle 1967, Hartle and Thorne 1968). Our preliminary results indicate that these 
improvements could bring a new information on the validity of EoS (Stuchlik et al. 
2007). The important role of the neutron star spin is demonstrated in the case of 
Circinus X-1 (Torok et al. 2010). 
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